Investigation of surface critical behavior of semi-infinite systems with cubic anisotropy. 
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The critical behavior at the special surface transition and crossover behavior from special to 
£Nj ' ordinary surface transition in semi-infinite n-component anisotropic cubic models are investigated by 

applying the field theoretic approach directly in d — 3 dimensions up to the two-loop approximation. 
The crossover behavior for random semi-infinite Ising-like system, which is the nontrivial particular 
case of the cubic model in the limit n — > 0, is also investigated. The numerical estimates of the 
resulting two-loop series expansions for the critical exponents of the special surface transition, surface 
qq . crossover critical exponent $ and the surface critical exponents of the layer, ati, and local specific 

' heats, an, are computed by means of Pade and Pade-Borel resummation techniques. For n < n c 

the system belongs to the universality class of the isotropic n-component model, while for n > n c 
the cubic fixed point is stable, where n c is the marginal spin dimensionality of the cubic model. 
i-C ' The obtained results indicate that the surface critical behavior of semi-infinite systems with cubic 

^ ' anisotropy is characterized by new set of surface critical exponents for n > n c . 

a: 

' I. INTRODUCTION 

a : 

^— > , 

The investigation of the critical behavior of real systems is a very important task of the condensed matter theory. 
The critical behavior in such systems as polymers, easy-axis ferromagnets, superconductors, as well as superfluid 

S 4 He, Heisenberg ferromagnets and quark-gluon plasma is described by isotropic 0(n) model with n = 0, 1, 2, 3 and 4, 
respectively, and is very well determined in the framework of different theoretical and numerical approaches. 

Investigation of the critical behavior of real cubic crystals is the subject of extensive theoretical work during past 
decades. In real crystals, due to their complex crystalline structure, some kind of anisotropy is always present. 
The simplest nontrivial crystalline anisotropy is a cubic one. A most typical model for description of the critical 
i behavior of such systems is the 0(n) model with a cubic interaction term ^ Y^7=i $ a dded to the usual % (| $ | 2 ) 2 ) 

termErmi This n-component cubic model is the particular case of an rrm-component model with cubic anisotropy at 
• m = 1. Anisotropic n-component cubic model exhibits different types of continuous and first-order phase transitions, 
depending on the number of spin components n, space dimensionality d, and sign of the cubic coupling constant vq. 
Cubic models are widely applied to the study of magnetic and structural phase transitions. In the limiting case of 
n — > 0, it describes the critical behavior of random Isinfelike systems. The case m = 1 and n — ► oo corresponds 
to the Ising model with equilibrium magnetic impuritiesu. Depending on the sign of the cubic coupling constant 
vq, two types of order are possible: along the diagonals of a hypercube in n dimensions (vq > 0) or along the easy 
CN| ! axes (vq < (X). In the latter case the system can undergo the first order phase transition, as was confirmed in some 
experimental In the present work we are concerned with the case vq > 0. 

The presence of surfaces, which are inevitable dnjreal systems, leads to additional complications. General reviews 
on surface critical phenomena are given in Refs.EHjij. A typical model to. study critical phenomena in real physical 
systems restricted by a single planar surface is the semi-infinite modelEl. There are different surface universality 
classes, defining the critical behavior in the vicinity of boundaries, at temperatures close to the bulk critical point 
(r = (T— T c )/T c — » 0). Each bulk universality class, in general, divides into several distinct surface universality classes. 
Three surface universality classes, called ordinary (cq — > oo), special (cq = c* p ) and extraordinary (co — > — oo), are 
O ■ relevant for our caseQutl; they meet at a multicritical point (m§,co) = (rap c ,c* ), w hich corresponds to the special 
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1.1) below, 
developed for pure isotropic 
and the systems with a ran- 



transition and is called the special point. The couplings mo and Co are defined in Eq.( 
TVipnry nf pritipal behavior of individual surface universality classes is very. -gel 
systemsffllijOlij'Lj, the systems with quenched surface-enhancement disorderu_j3tit 

dom quenched bulk disorder at the ordinary and the special surface transitions! 1 ^rl 18 !. General irrelevance-relevance 
criteria of the Harris type for the systems with quenched short-range correlated surface-bond disorder were predicted 
irJi.3 and confirmed by Monte-Carlo calculationsc3c3. Moreover, it was established that the surface critical behavior 
of semi-infinite systems with quenched bulk-disorder is characterized by the new set of surface critical exponents in 
comparison with the case of pure systemalJ'tS. 

Experimental systems are typically characterized by the parameters different from the fixed point values. In- 
vestigations of the crossover behavior between different surface universality classes, however,, hpte be e n - lim ited to 
pure, isotropic systems, in which quite rich and nontrivial crossover behavior has been found^llaEJO'EHlES. In the 
anisotropic systems the crossover behavior can be even more complex, because when the rotational symmetry is bro- 
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ken, the orientation of the surface may affect the surface critical phenomenaEil, and this effect can be different for 
different surface universality classes. 

In the present work we study the critical behavior at the special surface transition and crossover behavior from the 
special to the ordinary surface transitions, occuring in semi-infinite n-component model with cubic anisotropy. 

The effective Landau-Ginzburg- Wilson Hamiltonian of the model under consideration in the semi-infinite space is 
given by 



H{$) = J"dzj d d -\[\ | | 2 +\ m l | $ | 2 

+ ^oE^ + ^o(|0| 2 ) 2 )] + J d d -vi co 2 , (1.1) 



where 4>{x) is an n-vector field with the components <pi{x), i = 1, ...,n. Here m 2 , is the "bare mass", representing a 
linear measure of the temperature difference from the critical point value. The parameters uq and vq are the usual 



"bare" coupling constants uq > and vq > 0. In the case of replica limit n — > 0, the Hamiltonian (1.1) describes 
the critical, behav ior of semi-infinite Ising-like systems with a random quenched bulk disorder, when uq < and 
vq > (seeEZlE3£3). The constant cq is related to the surface enhancement, which measures the enhancement of the 
interactions at the surface. It should be mentioned that the rf-dimensional spatial integration is extended over a half- 
space IR^ = {x=(r, z) £ M d | r £ M d ~ 1 , z > 0}, bounded by a planar free surface at z = 0. The fields <fii(r, z) satisfy 
the Dirichlet boundary condition in the case of ordinary transition:- c£ujr, z) = at z = and th e N eumann boundary 
condition in the case of special transition: d z 4>i(v,z) — at z — 000. The model defined in ( [lTl| ) is translationally 
invariant in directions parallel to the external surface, z — 0. Thus, we shall use mixed representation, i.e. Fourier 
representation in d — 1 dimensions and real-space representation in z direction. 

The added cubic term breaks the 0(n) invariance of the model, leaving a discrete cubic symmetry. The model 



(1.1) has four fixed points: the trivial Gaussian, the Ising one in which n components are decoupled, the isotropic 
(O(n)-symmetric) and the cubic fixed points. The Gaussian and Ising fixed points are never stable for any number 
of components n. For isotropic systems, the 0(n)-symmetric fixed point is stable for n < n c , whereas for n > n c it 
becomes unstable. Here n c is the marginal spin dimensionality of the cubic model, at which the isotropic and cubic 
fixed points change stability, i.e. for n > n c , the cubic fixed point becomes stable. The 0(n)-symmetric fixed poit is 
tricritical. At n = n c , the two fixed points should coincide, and logarithmic corrections to the 0(n)-symmetric critical 
exponents are present. The calculation of the critical marginal spin dimensionality n c is the crucial point in studing 
the critical behavior in three-dimensional cubic crystals. Different results for n c have been published in a series of 
works in which different methods have been used. In the framework of the field-theoretical RG analysis the one-loop 
and three-loop approximations at e = 1 lead to the conclusion that n c should lie between 3 and 4E3i 2 J3, and the cubic 
ferromagnets are described by the Heisenberg model. On the other hand, by using the field theoretic approach directly 
in d = 3 dimensions u-a to the three-loop approximation, it has been found that n c — 2.9c3Ej. Similar conclusions 
were obtained in Ref.c3, where it was found that n c — 2.3. The calculations performed by Newman and RiedelEll 
with the help of the scaling-field method, developed by Goldner and Ridel for Wilson's exact momentum-space RG 
equations, h^e given for d = 3 the value n c = 3.4. Field-theoreijpaLjinalysis, based on the four-loop series in three 
dimensionstJtS and results of the five-loop e = 4 — d expansionooEj suggest that n c < 3. Recently a very precise 
six-loop result for the marginal|S»in dimensionality of the cubic model, n c = 2.89(4), was obtained in the framework of 
the 3D field-theoretic approaches. Thus, it was finally established that the critical behavior of the cubic ferromagnets 
is not described by the isotropic Heisenberg Hamiltonian, but by the cubic model, at the cubic fixed point. However, 
it was found that the difference between the values of the bulk critical exponents at the cubic and the isotropic fixed 
points is very small, and it is hard to measure this difference experimentally. The recently obtained results stimulate 
us to perform the investigation of the surface critical behavior of semi-infinite n-component anisotropic cubic model, 
and to determine corresponding surface critical exponents. 

The calculations are performed by applying the field theoretic approach directly in d = 3 dimensions, up to the 
two-loop order approximation. The numerical estimates of the resulting two-loop series expansions for the critical 
exponents of the special surface transition, and for the surface crossover exponent $ from the special to the ordinary 
transition arid surface criticaL-exponents of the layer, a±, and the local specific heat, an, are computed by means 
of the PadeEZI and Pade-BoreH resummation techniques for the cases n — 3, 4, 8 and for the case of n — > oo, which 
corresponds to the Ising model with equilibrium magnetic impurities. The crossover behavior from the special to 
the ordinary transition for random semi-infinite Ising-like system, which is the nontrivial particular case of the cubic 
model in the limit n — > 0, is also investigated. 
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II. RENORMALIZATION 

In order to investigate the critical behavior at the special surface transition in the semi-infinite n-component 
anisotropic cubic model and to calculate the crossover exponent <J>, we should consider correlation functions with 
insertions of the surface operator <jj?, 



N M Li 

i=i 3=1 



(2.1) 



1=1 



which involve N fields 0(xj) at distinct points Xj(l < i < N) off the surface, M fields <j>{tj,z = 0) = 4> s (i'j) at 

distinct surface points with parallel coordinates rj(l < j < M), and L\ insertions of the surface operator ^^(R;) 

(1 < ' < 'n). ... n 

The corresponding parallel Fourier transform of the full free propagator takes the forma 



G(p,z,z') = — 
2k 



Co 



C + K 



(2.2) 



with the standard notation kq = \/p 2 + m 2 ,. Here, p is the value of the parallel "momentum", i.e. the wave- vector, 
associated with d — 1 translationally invariant directions in the system. 

In the theory of semi-infinite systems the bulk field </>(x) and the surface field s (r) should be reparameterized 
by different uv-finite renormalization factorsSO, Z$(u,v) and Zl P (u,v). Thus, we have <fi — Z^ 2 (f)]i and <p s — 
Z^f 2 (Zl P ) 1 ^ 2 (/) St fj, where (j)R an d 4>s.b, ar( 3 the renormalized bulk and surface fields, respectively. Besides, introduction 

fiAn^l Gnr£an& r\-r\arct + r\r incprtinnc _, , , 

4>l = [ z <p 2 s 



of the additional surface operator insertions ^(j> 2 (R.i) requires additional specific renormalization factor Z$i, 



The corresponding renormalized correlation functions involving N bulk fields, M surface fields and L\ surface 
operators ^(p 2 (Ri) can be written as 



4 Wl) (0; rn, u, v, c) = Z -^ M y\z^)-^Z^G^ M ^ (0; m , u , v , c ). 



(2.3) 



In the present paper we concentrate our attention on the correlation function G(°' 2,1 )(0;m, u,v,c) involving two 
surface fields and a single surface operator insertion 0^(R;). 

The uv-singularities of the correlation function G^ N ' M,L ^ (0; p. u, v, c) can be absorbed through a mass shift mf, = 
m 2 + dm 2 and a surface-enhancement shift Co = c + 6c (seeH3). The renormalizations of the mass m, coupling 
constants u, fr-asd, the , renormalization factor Z^ are defined by standard normalization conditions of the infinite- 
volume theoryea'E9'EllE3. In order to absorb uv sin gulari ties located in the vicinity of the surface, a surface-enhancement 
s hift 8c is required (see Appendix 1). From Eq.( A1.3 ) and the expression for the renormalized correlation function 
(2.3) it follows that 



t^r 1 = z\ 



d[G^ 2 \Q-m 0: u^v^c Q )]- 1 



dco 



(2.4) 



cq— cq (c,m,w,u) 



It should be mentkmed that the renormalization factor Z^ = Z^Z^ is defined via the standard normalization 
condition (|A1.2|) (seeB@) 



d 



Z\l = ^-^[G^ 2 H P )Y 



= M 2»[GM(p)]-i. 

p^o p Op 



(2.5) 



It should be mentioned that all the Z factors in the d < 4 case have finite limits at A — ► oo (where A is the 
large-momentum cutoff) and depend on the dimensionless variables u and v. Besides, the surface renormalization 
factors Z sp and Zxi depend on both u,v and the dimensionless ratio c/m. The dependence on the ratio c/m plays a 
crucial role in the investigation of the crossover behavior from the special surface transition (c/m — ► 0) to the ordinary 
surface transition (c/m — > oo). 
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III. EXPANSION OF THE CORRELATION FUNCTIONS NEAR THE MULTICRITICAL POINT 

As was indicated before, the main goal of the present work is to investigate the critical behavior at the special 
surface transition and to perform the analysis of the scaling critical behavior between the special and the ordinary 
transition. In this connection let us consider the small deviation Aco = c — c* p from the multicritical point. The 
power expansion of the bare correlation functions G^ N ' M ' (0;mo,Uo,Vo, Co) in terms of this small deviation Acq has 
the form 



Li=0 



(3.1) 



Reexpressing the right-hand side of Eq.(p.l|) according to the Eq.(2.3) in terms of the renormalized correlation functions 
and renormalized variable Ac = [.2^,2 (it, v)]~ 1 Acq , we obtain 



z -iN +M)/ 2 {z s PrM/2GiN , M) (0 . ^ ^ Cq) 



(Ac) 



Li=0 



G^' M ' Li \0;m,u,v). 



(3.2) 



The above equation in a straightforward fashion defines the corresponding renormalized correlation functions, 
defined in the vicinity of the multicritical point 



G<f' M) (0; m, u, v, Ac) = ^-M/^M) (0; m ^ ^ VQjCQ y 

Ac/m 



(3.3) 



For dimensional reasons, we can introduce the dimcnsionlcss variable c — Aclm. Thus, the above correlation functions 
Gp M \0;m,u,v,,) satisfy the corresponding Callan-Symanzik equation^ 



d a , s d a , s 9 N + M . 

m dm + ^' V ' du + dv 2 V ' 



G^' M) (0;m, U ,«,c) = AG R , 



(3.4) 



where the inhomogeneous part AGr should be negligible in the critical region, similarly a s th is took place in the case 
of infinito^Si£S|tem field theory. The functions f3 u (u, v), (3 v (u, v) and 770 (it, v), appearing in (3.4), are the usual b ulk R jG. 
functionsL^lHcj. The resulting homogeneous equation differs from the standard bulk Callan-Symanzik equatiorc3c3'ED 
in that fashion that it involves the additional surface-related function r]i P and the term —[1 + r]c(u,v)]c-§=, where 



d 

77i p (u,u) = m— 



FP 



and 



r) E (u, v) = m 



(h 



FP 



dlnZ^iu, v) 
lnZ#i (u, v) = I3 U (it, v) g 2 - h /3 V (u, v) 



dv 

dlnZfyi (u, v) 
dv 



(3.5) 



(3.6) 



In the case Ac = we obtain the analog of the Callan-Symanzik equation for the correlation functions 
Gf^' M \0; m, u, v) at the special point cq = c* p . 

The symbol 'FP' indicates that the above value should be calculated at the corresponding infrared-stable fixed 
point (FP) of the underlying bulk theory. 



IV. ANALYSIS OF THE SCALING CRITICAL BEHAVIOR 



The asymptoti c sc aling critical behaviox-oLthe correlation functions can be obtained through detailed analysis of 
the CS equation (3.4), as was proposed irCj'cS and employed to the case of the semi-infinite systems ir£3E3o. In the 
critical region, at r = (T — T c )/T c — » 0, for the renormalization Z factors we obtain 



(4.1) 











Zfjj-2 
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where the variable identified with the inverse bulk correlation length £ _1 ~ t", as it is usually accepted in the 
massive field theorycJO'CJ. 

Substituting the last equation from (4.1) into the expressions for Ac and for the scaling variable c it is easy to 
obtain the following asymptotic forms 



Ac - m~^ u > v J Ac , 



Ac 



-u VE {u* ,v') 



(4.2) 



and 



-(i+i7 E («* ,v*)) AcQ) 



Ac , 



(4.3) 



where 



$ = v{\ +r)c(u*,v*)) (4.4) 

is the surface crossover critical exponent. The second equation in ( |4.3| ) explains the physical meaning of the surface 
crossover exp one nt as a va lue characterising the measure of deviation from the multicritical point. The second 
equations in (4.2) and ( [4.3] ) indicate non-analitic temperatu re d ependence of the renormalized surface-enhancement 
deviation Ac. Thus, from the analysis of the CS equation (3.4) we obtain the following asymptotic scaling form of 
the surface correlation function G^ 0,2 ' : 



G^ 2 \p- to , uo, v , c ) ~ m-^Gf^i^; 1, u* , v*, to^Acq) 

TO 

~r-^G(pt-*';l,T-*Atk,) ) (4.5) 
where 7^ = v(\ — 7711), is the local surface susceptebility exponent and 

V\\ P = VT + V4> ( 4 - 6 ) 

is the surface correlation exponent at the special surface transition. It is easy to see that the asymptotic scaling 
critical behavior of the surface correlation function for the semi-infinite ri-component systems with cubic anisotropy 
is characterized by the new crossover exponent <&(u*,v*), calculated at the cubic fixed point (u*,v*). 



THE FIXED-DIMENSION PERTURB ATI VE EXPANSION UP TO TWO-LOOPS. 



As was investigated previouslyBBIllil, the fixed-dimension </> 4 field-theoretic approacrS provides an accurate 
description of the surface critical behavior of semi-infinite systems. We apply this method to the analysis of the cubic 
anisotropic model (1.1) at the special surface transition and to the invctigation of the crossover behavior in the vicinity 
of the multicritical point Co = c* p . The surface correlation exponent 7/y P at the special transition can be obtained 

from Eq.(4.6), where if^ is defined by Eq.(3J5), and 770(14, v) is the standart bulk exponent 77^ = m-J^-lnZ '&\ „„. After 



\FP' 



performance of the mass- and surface-enhancement renormalization and carrying out the integration of the corre- 
sponding Feynman integrals with subsequent execution of the standart vertex renormalizations of bare dimensionless 



parameters uq =r^(l + 6^ + v), vq = v(l + |w + 2u) (where uq — 7*0/877777 and vq = vo/87tto), by analogy to that 
as it was done mt3 for random semi-infinite model, we obtain in the case of the cubic anisotropic model the following 
expression for the surface correlation exponent at the special transition 



v7( u , v ) 



n+2 v 
~2(n+8) U ~ 6 



(5.1) 



(n+8) 2 V ' 9 K ' n+8 



A(n)uv. 



Here A(n) is defined as 



^ „ . n—10 n + 2,, 9 

A(n) = 2A H 1 77 2 2 - ln2), 

v > 48 6 v ; ' 



(5.2) 



and the renormalized coupling constants u and v, are normalized in the standard fashion, 77=2^77 and v=%v. 
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According to the Eq.(4.4) and Eq.([3.6]), the calc ulati on of the crossover critical exponent <!> is connected with the 
calculation of the rcnormalization factor Z^ via (2.4). With that end in view we can rewrite the normalization 
condition (A. 1.1) in the form 



Z\\ [G (0,2) (0; to , "o, «o, c )] 1 = to + c 



(5.3) 



for the inverse unrenormalized surface correlation function [G^°' 2 ^(0)] _1 . It is easy to see that b y di fferentiation of 
the above mentioned normalization condition with respect to cq and by taking into account Eq.(2.4) we obtain for 
the renormalization factor Z^i the equation 



where cq 



5c and 



dcp 
dc 



5c = (Zn 1 — l)(m + c) + (To(0; to, c = c + 5c). 



(5.4) 



(5.5) 



Here-jge/O; to, cq) — o\ + 02 + 03 + 04 denotes the sum of loop diagrams of all orders in [G^ 0,2 ^(0; to, uq, Vo, cq)] 1 
(seelSES). Among them u\ corresponds to the one-loop graph, 02 denotes the melon-like two-loop diagrams 



02 



G 
G 
G 



2k 



to 2 d 



(5.6) 



k 2 =0 



(T3 and (7 4 represent the one-particle reducible and one-particle irreducible two-loop diagrams in 



[G(°' 2 )(0; 



metry factors) 



to, Mo, ^o, Co)] , respectively. These graphs have their correponding weights (arise from the standard sym- 



h_ with $) = n + 
2 1 3 



(0) 



with 4 0) = - 



■ Uq + V . 



ul + vl + 2v u , 



,(0) 



f (0) 



^ and W i th 4°) = tf = 



(5.7a) 
(5.7b) 
(5.7c) 



The full lines with notes "G" denote the full free propagator (2.2) a nd s igns "-" denote t he fr ee bulk propagators, 
which are associated with the first term in the full free propagator (2^2). The equation (5.4) can be resolved by 
using the method of sequential iteration. As the result, at the second order of the perturbation theory, we obtain the 
renormalization factor Zm2 in terms of the new renormalized coupling constants u and v, 



Zj,2(u,v) = 1 + 



where C(n) is a function defined as 



> 1 1 n + 

'-{ln2 - -)u + (ln2 - -)v H — 

v 4 y v A' 3 



-C(n)u 2 + 2C(n)uv + G(l)^ 2 , 



C{n) =A-B- -ln2 



llnH- 2n+1 



12 



(5. 



(5.9) 



and A = 0.202428, B = 0.678061 are integrals originating from the two-loop melon-like diagrams. Combining 
the renormalization factor Z^ with the one-loop pieces of the j3 functions 0n(u, v) — —u(l — [(n + 8)/6]u — v) and 



(5.10) 



fiv(u, v) = —v(l — iv — 2-u), according to Eq.(p.6[), we obtain the desired series expansion for r/g, 



7] S (U,V) 



-2^±|Gn2 - \)u - \{ln2 - \)v - 8(3^^ 
n + 8 4 3 4 (n + 8) z 



D(n)u 2 



2D{n) D(l) 



where 



D(n) 



A-B 



n + 2, ,„ n + 1, „ 17n + 22 

ln 2 2 hi2 H . 

3 2 96 



(5.11) 



Eq.(p^|) gives a result for the model with the effective Hamiltonian of the Landau- Ginzburg- Wilson type with cubic 
anisotropy in the semi-infinite space ( ji . 1| ) , with a general number n of order parameter components. 
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The knowledge of ry E gives an access to the calculation of the crossover critical exponent <f> via the scaling relation 



(4.4). Additionally, we can calculate the critical exponents ot\ and an of the layer and specific heats via the usual 
scaling relationsQ 

ai = a + i/-l + $ = l-i/(d-2-7?c), a n = a + v - 2 + 2$ = -v\d - 3 - 2r) s ]. (5.12) 

The above critical exponents should be calculated for different n (n = 3, 4, 8, and n oo) at the standard infrared- 
stable cubic fixed (FP) points of the underlying bulk theory, as it is usually accepted in the massive field theory. 
As it was mentioned above, in the cases n < n c the cubic ferromagnets are described by the Heisenberg isotropic 
Hamiltonian at the 0(n)-symmetric fixed point. 



In the case of the replica limit n — > we can obtain from (5.10) the series of 7/g and corresponding expressions 
for the surface crossover critical exponent. $ r , and critical exponents a\ and a r xl of the layer and specific heats for 
semi-infinite random Ising-like systems (c3) 



VI. NUMERICAL RESULTS 



In order to obtain the full set of surface critical exponents at the special transition and to calculate the surface 
crossover e xpon ent $ from the special to the ordinary transition in systems with cubic anisotropy we substitute the 
expansion ( |5.l| ) for rj^ p into the standart scaling-law expressions for the surface exponents (see Appendix 2) and the 
expansion ( |5.10| ) for r/ e into the scaling relations (4.4) and (5.12). 

For each of the above mentioned surface critical exponents of the special transition and for the crossover exponent 
$ we obtain |t al d = 3 a double series expansion in powers of u and v, truncated at the second order. As it is 
knowr£3H'L3E3, power series expansions of this kind are generally divergent due to a nearly factorial growth of 
expansion coefficients at large orders of perturbation theory. In order to perform the analysis of these perturbative 
series expansions and to obtain accurate estimates of the surface critical exponents—a powerful resummation procedure 
must be used. One of the simplest ways is to perform the double Pade-analysi£3. This should work well when the 
series behave in. lowest orders " in a convergent fashion" . Another way is to perform for these series the double Pade- 
Borel analysisE3. jThe Pade-Borel resummation procedures are posible to use only in the case when the series are 
alternating in signo. 

The results of our calculations of the surface critical exponents of the special transition for various values of 
n = 3,4, 8, oo at corresponding cubic fixed points are presented in Tables 1-5. Unfortunately, the second (p = 2) 
order analysis of perturbative serieaEj gives the cubic fixed point with coordinates uo — 1.5347 and vq = —0.0674 at 
n — 3 for the 3D model. This corresponds to ordering along the easy axes, because (vq < 0). The analysis of the 
eigenvalues of the stability matrix shows that in the frames of the two-loop approximation the cubic fixed point at 
n = 3 is unstable and the 0(n)-symetric fixed po int is stable. BytU-the estimations of the marginal-stm .dimensionality 
of the-eubic model n c in the frames of three-loopc§Ej, four-loop£3E3, five-loop e = 4 — d expansioiSE-J'Ea and six-loop 
studyE3 show that the cubic ferromagnets are not described by the Heisenberg isotropic Haiailtonian, but by the 
cubic model at the stable cubic fixed point. Higher precision six-loop field-theoretical analysisE3 give the value of the 
marginal spin dimensionality of the cubic model equal to n c — 2.89(4). In accordance with this we used the cubic 
fixed point of the higher p — 3 order of perturbative series for obtaining the set of surface critical exponents at n — 3. 
For estimation of the reliability of the obtained results we performed calculations at the cubic fixed point of the p = 6 
order in Table 2. The obtained results indicate that the difference in the ways of the (3 functions resummation have 
no essential influence on the values of the surface critical exponents and that the results obtained in the frames of the 
two-loop approximation are stable and reliable. The surface critical exponents of the special transition for n = 4, 8 
and n — > oo were calculated at the standard infrared-stable cubic fixed (FP) points of the underlying bulk theory, as 
it is usually accepted in the massive field theory. 

The Tables 6-9 represent the surface crossover exponent $ and the surface critical exponents of the layer a\ and 
local an specific heats for different values of n — 3, 4, 8, oo at the corresponding cubic fixed points. 

The quantities 0\jOi and Ou/02i represent the ratios of magnitudes of first-order and second-order perturbative 
corrections appearing in direct and inverse series expansions. The larger (absolute) value of these ratios indicate the 
better apparent convergence of truncated series. 

The values [p/q] (where p, q = 0,1) are simply Pade approximants which represent the partial sums of the direct and 
inverse series expansions up to the first and the second order. The nearly diagonal two- variable rational approximants 
of the types [11/1] and [1/11] give at u — or v — the usual [1/1] Pade approxirnantO. The results of Pade-Borel 
analysis of the direct R and the inverse series expansions give numerical estimates of the surface critical exponents 
with a high degree of reliability. As it is easy to see from Tables 1-5, the most reliable estimate is obtained from the 
inverse series expansion for the surface critical exponent rj±, which represent the best convergence properties. Thus 
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we obtain for r]±: —0.121 at n = 3; —0.139 at n — 4; —0.143. at n — 8; —0.141 at tw^ oo. Substituting these values 
of r)± together with the standard bulk values^ (v = 0.6970 (or 0.706(6X.at p = QB) (n = 3); 0.711 (n = 4); 0.717 
(n = 8); 0.715 (*h-> oo))E3 and 77 (77 = 0.024E2 (or 0.0333(26) at p = 6E2) (n = 3); 0.026 (n = 4) ; 0.025 (n = 8); 
0.025 (n — > oo))E3 at <i = 3 into the surface scaling relations (see Appendix 2) we have obtained the remaining 
critical exponents that are present in the last columns of Tables 1-5. The deviations of these estimates from the other 
estimates of the table give a rough measure of the achieved numerical accuracy. 

The double Pade and Pade-Borel analysis of the series for r] s , a±, an and presented in Tables 6-9, was performed 
in a similar way. In this case the situation is more complicated, because the series of these critical exponents exhibit 
bad convergence properties. A very similar situation took place in the analysis of perturbation expansions for the 
surface critical exponents in isotropic systemstj. As it is easy to see, we obtain the most reliable estimates for the 
surface critical exponent an- The result of substituting of this value an (an = — 0.33 1 (n = 3); —0.374 (n = 4); 
—0.384 (n = 8); —0.380 (n — » 00)), together with the value of v, into scaling laws (4.4) and ( 5.1 2| ) is presented in 
the last columns of Tables 6-9. For evaluation of the reliability of the obtained results we have performed additional 
calculation ofj-the surface critical exponents on the basis of these an and the corresponding six-loop perturbation 
theory resultsEj for bulk critical exponent of the correlation length v = 0.706(6) (n = 3); 0.714(8) (n = 4); 0.712(6) 
(n = 8); 0.708(8) (n — > 00). The results of calculation are presented in Table 10. As it is easy to see, these results 
differ very little from the results presented in the last columns of Tables 6-9. It indicates good stability of the results 
obtained in the frames of the two-loop approximation scheme. 

The obtained results for surface critical exponents of semi-infinite model with cubic anisotropy, calculated at the 
cubic fixed. poiat-arc different from the results for surface critical exponents of standart semi-infinite n-componemt 
model (seeaolij'Lj). For example, the difference in the case of d = 3 at n = 3 is about 7.5% for m and 6.6% for rj^E3. 



If n < n c , the cubic fixed point is unstable and the cubic term in the Hamiltonian (1.1) becomes irrelevant. In this 
case the isotropic fixed point is stable and the system is described by the simple O (n)-jsvmmetric model in 3D. The 
corresponding surface critical exponents can be calculated from the seri es, p resented inE3. 

As was indicated previously, in the limit n — > 0, the cubic model (1.1) with uo < and vo > describes the 
semi-infinite Ising-like systems with, random bulk disorder. The investigation of the surface special transition for such 
kind of systems was presented inlla. For the value of rj^, for the surface crossover exponent <fr r and for the surface 
exponents of the layer a\ and locai-a^ specific heat, in the limit n — > at th&,standard infrared-stable random fixed 
point u* = -0.60509 v* = 2.39631E3 of the underlying bulk theory we obtainEl 



4 = -0.164, 
a r n = -0.222 



a\ = 0.211, 
$ r = 0.567. 



(6.1) 



These values are different from the surface critical exponents of the pure semi-infinite Ising-like systems and indicate 
that in a system with random bulk disorder the planar boundary is characterized by a new set of the surface critical 
exponents. 



VII. SUMMARY 



We have studied special surface transition and crosover critical behavior in the vicinity of the multicritical point for 
3D semi-infinite systems with cubic anisotropy by applying the field theoretic approach directly in d — 3 dimensions 
up to the two-loop approximation. We have performed a double Pade and Pade-Borel analysis of the resulting 
perturbation series for the surface critical exponents of the special transition, the surface crossover exponent $ and 
the critical exponents of the layer ai and local specific heats an for various n — 3,4,8, 00, in order to find the best 
numerical estimates. 

We find that at n > n c the surface critical exponents of the special transition in semi-infinite systems with cubic 
anisotropy belong to the cubic universality class and the asymptotic scaling critical behavior of the surface correlation 
function in the vicinity of the multicritical point is characterized by the new crossover exponent of the cubic universality 
class <&(u*, v*). 

Besides, we present the results of investigation of the crossover critical behavior between special and ordinary 
transition in random semi-infinite systems, by taking the limit n — > 0£J. In this case we calculate the surface critical 
exponents a\, and the surface crossover exponent <i> r and confirm that they belong to the universality class of the 
random model. 

The further theoretical investigation of the asymptotic surface critical behavior of semi-infinite cubic systems will 
be highly desirable in the framework of higher-order (three-loop) RG approximations. In particular, it would be very 
interesting to study the case of the 3D cubic crystal with n = 3. 

We suggest that the obtained results could stimulate further experimental and numerical investigations of the 
surface critical behavior of random systems and systems with cubic anisotropy. 
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Appendix 1 

We require that the following surface normalization conditions take placeHH : 

Gf 2 \0;m,u,v,c) = ^—, (Al.l) 
n m + c 



dGft ' 2 ^ (p; m, u, v, c) 



dp 2 



p=0 



2m(m + c) 



2 ' 



(A1.2) 



G^ 2 ' 1) (0;m, M , U , C ) = ^-l-^ (A1.3) 

where the correlation function G*- ' 2 ' 1 -' contains the insertion of the surface operator \<ft 2 s - 

The Eq.( A1.5| ) is motivated by the fact that the bare correlation function G^ ' 2,1 ^(0; m 2 ,, uq, vq, cq) may be written 
as a derivative — ^-G^ 0,2 ) (0; m§, uq, vq, cq). This equation simplifies considerably the calculation of the correlation 
function with insertions of the surface operator i 



2 



It is easy to see from (Al.l) that the special point is located at m = c = 0, because at this point the divergence of the 



bulk and the surface correlation length and susceptibility is observed. At c = the surface normalization conditions 
are simplified and yield cq = c* p . This point corresponds to the multicritical point {itiq c , c* p ) at which special transition 
takes place. On the other hand, the above mentioned equation implies also that the surface correlation length and 
the susceptibility are finite at the ordinary transition, because in this case we have c > when m — > 0. This latter 
case corresponds to the situation when the surface remains " noncritical" at the bulk transition temperature. 



Appendix 2 



The individual RG series expansions for other critical exponents can be derived through standard surface scaling 
relational] with d = 3 



T) + 1 


)\\ 


2 




V -{d- 


-2 + i7||), 


--y{l- 


-*?ll)> 


i/(2- 




& 




A 


d + 2-T] 


Ji ~ 


d-2+m' 


Ai 


d- m 


0i ' 


d~2 + r]\ 



Each of these critical exponents characterizes certain properties of the cubic anisotropic system near the external 
surface, in the vicinity of the critical point. The values v, rj, and A = v(d + 2 — 77) /2 are the standard bulk exponents. 
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TABLE I: Surface critical exponents of the special transition for d — 3 up to two-loop order at the cubic fixed point (of order 
p = 3) u* = 1.348, v* = 0.074 at n = 3. 



exp 


L/ 1 

o ? 


W 1 7. 

C 7 ,i 


[0/0] 


[1/0] 


[o/i] 


[2/0] 


[0/2] 


[H/1] 


[1/H] 


R 


R 




V\\ 


-4.0 


14.7 


0.00 


-0.319 


-0.242 


-0.239 


-0.254 


-0.255 


-0.255 


-0.260 


— 


-0.266 


Ai 


10.8 


-4.4 


0.75 


1.069 


1.218 


1.098 


1.078 


1.101 


1.102 




1.109 


1.138 


V-L 


-3.1 


-6.0 


0.00 


-0.159 


-0.137 


-0.107 


-0.117 


-0.121 


-0.120 


-0.124 


-0.121 


-0.121 


Pl 


0.0 


0.0 


0.25 


0.25 


0.25 


0.253 


0.253 










0.256 


7n 


11.8 


-4.3 


0.50 


0.819 


0.968 


0.846 


0.823 


0.848 


0.849 




0.857 


0.882 


7i 


12.6 


-3.1 


1.00 


1.398 


1.662 


1.430 


1.372 


1.433 


1.433 




1.450 


1.478 


Si 


-3.8 


-1.7 


5.00 


6.593 


7.339 


6.236 


5.763 


6.262 


6.266 


6.287 


6.359 


6.779 


5ii 


-5.8 


-1.7 


3.00 


4.275 


5.217 


4.057 


3.622 


4.090 


4.095 


4.101 


4.186 


4.450 



TABLE II: Surface critical exponents of the special transition for d = 3 up to two-loop order at the cubic fixed point (of order 
p = 6) u* = 1.321(18), v* = 0.096(20) at n = 3. 



exp 


Oi 

o 2 


Oii 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[H/1] 


[1/H] 


R 


R' 1 






-4.0 


14.9 


0.00 


-0.316 


-0.240 


-0.237 


-0.252 


-0.254 


-0.253 


-0.258 




-0.275 




10.9 


-4.5 


0.75 


1.066 


1.212 


1.095 


1.075 


1.098 


1.099 




1.106 


1.156 


V-L 


-3.1 


-6.0 


0.00 


-0.158 


-0.137 


-0.107 


-0.117 


-0.120 


-0.120 


-0.123 


-0.121 


-0.121 


Pi 


0.0 


0.0 


0.25 


0.25 


0.25 


0.253 


0.253 










0.256 


7n 


11.9 


-4.3 


0.50 


0.816 


0.962 


0.843 


0.821 


0.845 


0.846 




0.854 


0.900 


7i 


12.7 


-3.2 


1.00 


1.395 


1.654 


1.426 


1.370 


1.429 


1.430 




1.446 


1.497 


Si 


-3.8 


-1.7 


5.00 


6.581 


7.312 


6.248 


5.765 


6.255 


6.260 


6.279 


6.351 


6.853 


5u 


-5.9 


-1.7 


3.00 


4.265 


5.187 


4.050 


3.624 


4.083 


4.089 


4.094 


4.178 


4.520 



TABLE III: Surface critical exponents of the special transition for d — 3 up to two-loop order at the cubic fixed point (of order 
p = 2) u* = 1.064, v* = 0.520 at n = 4. 



exp 


Oi 

o? 


0?A 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[H/1] 


[1/H] 


R 


R- 1 


f(v±,^ v) 


■n\\ 


-4.6 


7.3 


0.00 


-0.353 


-0.261 


-0.277 


-0.286 


-0.294 


-0.288 


-0.299 




-0.304 


Ai 


7.2 


-4.7 


0.75 


1.103 


1.295 


1.152 


1.134 


1.157 


1.165 




1.173 


1.175 


V-L 


-3.5 


-8.8 


0.00 


-0.176 


-0.150 


-0.125 


-0.135 


-0.139 


-0.138 


-0.143 


-0.139 


-0.139 


01 


0.0 


0.0 


0.25 


0.25 


0.25 


0.254 


0.254 










0.247 


7n 


7.2 


-4.7 


0.50 


0.853 


1.045 


0.902 


0.883 


0.906 


0.915 




0.924 


0.927 


7i 


8.0 


-3.2 


1.00 


1.441 


1.789 


1.496 


1.432 


1.501 


1.512 




1.532 


1.521 


*i 


-4.3 


-1.7 


5.00 


6.764 


7.725 


6.762 


5.861 


6.454 


6.487 


6.479 


6.596 


7.147 


Sii 


-7.8 


-1.7 


3.00 


4.411 


5.664 


4.231 


3.699 


4.262 


4.303 


4.271 


4.413 


4.747 
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TABLE IV: Surface critical exponents of the special transition for d = 3 up to two-loop order at the cubic fixed point (of order 
p = 2) u* = 0.525, v* = 1.146 at n = 8. 



exp 


Oj 

o 2 


Oii 

o M 


[0/0] 


[1/0] 


[o/i] 


[2/0] 


[0/2] 


[11/1] 


[1/H] 


R 


R- 1 


f(r)±, v, rj) 


'''II 




fi 9 
u . z 


n no 


-yJ.OOO 


9fi9 


-U. ZrOU 


9Q9 


-U.OU4 


9Q9 


-u.ouo 




n ^1 1 


Ai 


6.5 


-5.0 


0.75 


1.105 


1.300 


1.160 


1.146 


1.164 


1.178 




1.187 


1.187 


VX 


-3.8 


-11.6 


0.00 


-0.177 


-0.151 


-0.131 


-0.140 


-0.145 


-0.142 


-0.148 


-0.143 


-0.143 


Pi 


0.0 


0.0 


0.25 


0.25 


0.25 


0.257 


0.257 










0.247 


7n 


6.4 


-5.1 


0.50 


0.855 


1.050 


0.911 


0.898 


0.913 


0.931 




0.940 


0.940 


7i 


7.1 


-3.3 


1.00 


1.444 


1.798 


1.506 


1.447 


1.509 


1.532 




1.552 


1.537 


Si 


-4.8 


-1.8 


5.00 


6.775 


7.751 


7.125 


5.920 


6.505 


6.570 


6.528 


6.676 


7.221 


Su 


-9.6 


-1.7 


3.00 


4.420 


5.695 


4.272 


3.750 


4.302 


4.382 


4.310 


4.490 


4.806 



TABLE V: Surface critical exponents of the special transition for d = 3 up to two-loop order at the cubic fixed point (of order 
p = 2) u* = 0.201, v* = 1.508 at n -» oo. 



exp 


Oi 
Oj 


o„ 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[H/1] 


[1/H] 


R 


R- 1 


f(71x,v, V) 


m 


-5.0 


6.5 


0.00 


-0.352 


-0.260 


-0.282 


-0.289 


-0.301 


-0.288 


-0.306 




-0.307 




6.7 


-4.9 


0.75 


1.102 


1.293 


1.154 


1.140 


1.157 


1.173 




1.182 


1.182 


TlX 


-3.7 


-10.7 


0.00 


-0.176 


-0.150 


-0.129 


-0.138 


-0.144 


-0.140 


-0.147 


-0.141 


-0.141 


Pi 


0.0 


0.0 


0.25 


0.25 


0.25 


0.258 


0.258 










0.248 


711 


6.6 


-5.0 


0.50 


0.852 


1.043 


0.905 


0.891 


0.906 


0.926 




0.935 


0.935 


7i 


7.4 


-3.3 


1.00 


1.440 


1.785 


1.499 


1.440 


1.501 


1.526 




1.546 


1.531 


Si 


-4.7 


-1.8 


5.00 


6.759 


7.714 


7.204 


5.906 


6.492 


6.564 


6.515 


6.669 


7.179 


Su 


-9.1 


-1.7 


3.00 


4.407 


5.651 


4.252 


3.738 


4.287 


4.376 


4.296 


4.483 


4.772 
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TABLE VI: Surface critical exponents involving the RG function r\c for the case n = 3, d = 3 at the cubic fixed point (of order 
p = 3) u* = 1.348, w* = 0.074. 



exp % [0/0] [1/0] [0/1] [2/0] [0/2] [11/1] [1/11] [R] R^ f(au,u, v ) 
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TABLE VII: Surface critical exponents involving the RG function r/c for the case n = 4, d — 3 at the cubic fixed point (of order 
p = 2) u* = 1.064, v* = 0.520. 
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